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Following the recent success of anisotropic hydrodynamics we propose a new, general prescription 
for the hydrodynamics expansion around an anisotropic background. The anisotropic distribution is 
fixing exactly the complete energy-momentum tensor, just like the effective temperature is fixing the 
proper energy density in the ordinary expansion around local equilibrium. This means that momen¬ 
tum anisotropies are already included at the leading order, allowing for large pressure anisotropies 
without the need of a next to leading order treatment. The first moment of the Boltzmann equation 
(local four-momentum conservation) provides the time evolution of the proper energy density and 
the four velocity. Differently from previous prescriptions, the dynamic equations for the pressure 
corrections are not derived from the zeroth or second moment of the Boltzmann equation, but they 
are taken directly from the exact evolution given by the Boltzmann equation. We check the effec¬ 
tiveness of this new approach by matching with the exact solution of the Boltzmann equation in 
the Bjorken limit with the collisional kernel treated in relaxation time approximation, finding an 
unprecedented agreement. 

PACS numbers: 12.38.Mh, 24.10.Nz, 25.75.-q, 51.10.+y, 52.27.Ny 

Keywords: relativistic heavy-ion collisions, quark-gluon plasma, anisotropic dynamics, viscous hydrodynam¬ 
ics, Boltzmann equation, relaxation time approximation, RHIC, LHC 


I. INTRODUCTION 

Relativistic viscous hydrodynamics has been an important tool for understanding the evolution of matter in the 
extreme conditions of heavy-ion collisions at RHIC (Relativistic Heavy-Ion Collider) and the LHC (Large Hadron 
Collider). Therefore a large interest has been triggered in the development of the hydrodynamic framework [1-27]. 
Despite the succesfull use of viscous hydrodynamics, some theoretical shortcoming have to be properly addressed. 
Several studies pointed out that viscous corrections combined with rapid longitudinal expansion (caused by large 
longitudinal gradients) result in a substantial pressure asymmetry at early times, even at a very low viscosity to 
entropy density ratio, 4ttt)/S ~ 1. Similarly, a variety of microscopic models (string models, color glass condensate, 
pQCD kinetic calculations) predict large momentum anisotropies at early times. On the other hand, in the limit of 
infinitely strong coupling, where the AdS/CFT correspondence can be used as an effective model, one finds a similar 
situation — a large difference between the pressures in the transverse plane and in the longitudinal direction [28, 29], 
persisting for a significant fraction of the evolution. Large pressure anisotropies are a cause for concern, since viscous 
hydrodynamics treats pressure corrections in a perturbative manner and, indeed, the presence of very large shear 
corrections (of the order of the isotropic pressure) may lead to unphysical results such as the negative longitudinal 
pressure and/or negative one-particle distribution functions. 

A new approach to treat these problems is anisotropic hydrodynamics (aHydro) [30-42]. The main feature of this 
approach, in contrast with the conventional approach, is to treat the expected high degree of pressure anisotropy 
of the produced matter already in the leading order of the hydrodynamic expansion. In the original formulation of 
aHydro, the leading order includes a single anisotropy parameter taking into account a large difference between the 
longitudinal and transverse pressures. This approach is unable to reproduce the pressure anisotropy in the transverse 
plane, which is generated by the radial flow [40]. However, non-trivial transverse dynamics is essential for the correct 
description of transverse collective behavior, like the elliptic flow, i> 2 , suppression. 

In Ref. [43], Bazow, Heinz, and Strickland extended the aHydro framework to the (2+l)-dimensional case, which is 
longitudinally boost invariant but non-homogeneous in the transverse plane. This approach has been recently extended 
to massive (non-conformal) systems in Ref. [44], In Refs. [43, 44], however, the transverse pressure asymmetries (and 
part of the longitudinal pressure correction) are included only through next to leading order corrections. On the 
other hand, it has been demonstrated (see Ref. [45-48]) that it is possible to generalize the background to the (1+1)- 
dimensional case (boost invariant and cylindrically symmetric radial flow). Very recently, we have proposed a set of 
equations for anisotropic hydrodynamics that is consistent with the full (3+l)-dimensional expansion [49]. 
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In all of the previous works on anisotropic hydrodynamics there was an ambiguity in the choice of the dynamical 
equations of motion, both for the leading order treatment and the next to leading order of the hydrodynamics 
expansion. Besides the local four-momentum conservation (first moment of the Boltzmann equation), at least some of 
the remaining equations have been derived from the zeroth moment (typically describing particle creation), the second 
moment of the Boltzmann equation or both. However, in principle any other moment of the Boltzmann equation could 
provide the necessary equations. This is not a new problem, since in the ordinary viscous expansion there is exactly 
the same ambuguity [50]. As pointed out in Ref. [50], this problem has already been solved for the hydrodynamics 
expansion around local equilibrium. The Boltzmann equation can be used to provide directly the evolution equation 
for the moments of the Boltzmann distribution (not to confound with the moments of the Botzmann equation). Among 
them, for instance, there are the shear pressure corrections tt ^ v , and the higher moments which appear in the exact 
evolution of n itself. This approach has been very successful for second order viscous hydrodynamics [51, 52], 
however, it can not be extended directly to the anisotropic expansion. Since the anisotropic background has at least 
one extra free degree of freedom (anisotropy parameter, contributing to the shear viscous corrections), the exact 
evolution from the Boltzmann equation (for instance the evolution of will not be able to fix both the leading 
order and the next to leading order contribution (in previous works denoted as + 7f M ", with 7rdepending 

only on the anisotropic background [43, 44, 49]). 

In this paper we propose to use a generalized version of the Landau matching in order to determine extra degrees 
of freedom related to anisotropy. Just like the effective temperature fixes the proper energy density at all orders in 
the viscous hydrodynamics expansion, we propose that the extra parameters in the anisotropic distribution will fix 
the pressure corrections in the anisotropic expansion (in this case n 1 * 1 ' = tt)) 1 ' at all orders by definition). In order to 
illustrate our new method, we consider a particular application of the proposed prescription, namely, we deal with the 
leading order of the expansion with the collisional kernel treated in the relaxation-time approximation, and we test it 
against the exact solutions in the Bjorken flow limit. The agreement is significantly improved compared both to the 
best prescription for leading order anisotropic hydrodynamics, and to the best prescription for second order viscous 
hydrodynamics (which is proved to be the close-to-equilibrium limit of this novel prescription), without the need to 
consider the next to leading order. 

The paper is organized as follows: In the next Section we introduce the hydrodynamics expansion and discuss the 
anisotropic case, including the generalized anisotropic matching procedure. The selection of the dynamic equations 
is discussed in Sec. III. Section IV contains the analysis of the specific case of the collisional kernel treated in the 
relaxation-time approximation, the close-to-equilibrium limit, and the numerical results. Summary and conclusions 
are presented in the last Section. Throughout the paper the natural units are in use, where c = h = ks = 1, and 
the metric tensor has the signature (+, —, —, —). The contraction of vectors (and rank two tensors) can be expressed 
both by repeated indices in the upper and lower position U^T^Uv and by a dot U -T -U. Details of the calculations 
for the 0 + 1 dimensional case are presented in the Appendix. 


II. HYDRODYNAMICS EXPANSION 

The purpose of this work is to extract hydrodynamics from the relativistic Boltzmann equation. For mathematical 
simplicity we consider only a single fluid of massive particles, without external fields, in flat space-time. Therefore, 
the Boltzmann equation reads 


P fl d ll f{x,p) = -C[f], (1) 

where f(x,p) is the particle distribution and C[f] is the collisional kernel, in the most general case a functional of the 
distribution /. Using the kinematic definition of the energy-momentum tensor we can write the local four-momentum 
conservation as the first moment of the Boltzmann equation (1). Indeed, with the energy-momentum tensor T^ 
defined as 


T^(x) = J dPp^p" f(x,p), 

the first moment reduces to the ordinary divergence-free condition on the energy-momentum tensor 


( 2 ) 


0 = J dPp u C[f] = J dPp v p ■ Of = d t , 


dPp> J p IJ f = d tl T^. 


( 3 ) 


3 


We understand that f dP in the last formulas is the Lorentz invariant momentum measure J dP = f d 3 p/E p , with 
E p = y/m 2 + p 2 the on-shell energy of particles with mass m. 

Equation (3) is the fundamental equation of hydrodynamics, governing macroscopic evolution of energy and mo¬ 
mentum densities 1 . Unfortunately (3) contains only four independent equations out of the ten degrees of freedom of 
the energy-momentum tensor, therefore it cannot be solved directly. 


A. Viscous hydrodynamics expansion 

Traditionally, in order to extract viscous hydrodynamics from an underlying kinetic theory, that is, in order to 
find the necessary extra equations needed to close the system, one expands the distribution function around the local 
equilibrium 


f{x,p) = f eq .{x,p) + Sf(x,p). 


For Boltzmann-like statistics, the local equilibrium distribution reads 


fe q .(x,p) 


k exp 


PyU^(x) 

T(x) 


( 4 ) 

( 5 ) 


which is the point dependent form of the global equilibrium distribution. It is assumed at this stage that the deviation 
from equilibrium is small and it can be treated in a perturbative manner. For a detailed discussion about the meaning 
of a perturbative treatment of Sf, and about different approaches to obtain these extra equations see Ref. [50]. 

In the formula (5) we have not yet defined the fluid four-velocity I/ M and the effective temperature T. A common 
definition for U M is the Landau definition 2 . The four-velocity is the the time-like eigenvector of the energy-momentum 
tensor 


UpT^ = £ U v , (6) 

where £ is the proper energy density, namely the energy density seen from the local rest frame. In other words, the 
requirement in Eq. (6) is equivalent to the condition on Sf 

AT a0 = J dP(p- U)pM Sf = 0. (7) 

With the notation A "^’" = A{) A — v — we are indicating the projection orthogonal to the four-velocity for the index 
/i, with A ^ being the projector 


A^ = g^ - U^UA (8) 

The usual definition of the effective temperature is that it gives the correct proper energy density 

£ = U il U v = J dP (p- U) 2 f = J dP(p ■ f/) 2 /eq., (9) 

which is generally referred to as the Landau matching. This is implicitly an additional constraint on the possible 
forms of Sf , before making any approximations, namely 

J dP (p ■ U) 2 Sf = 0. (10) 

We can present at this point the most general decomposition of the stress-energy tensor consistent with the Landau 


1 In presence of additional conserved charges like the baryon number, one should add the relative equations, but in the present paper we 
consider the simplest case. 

2 Another popular definition is the Eckart one. It must be noted that these are not the only two possible definitions (see, for instance, 

Refs. [54—56] for extended discussion). 
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prescription for the four-velocity 


= £ U^U V - (p eq . +n) + 7 T^. (11) 

In this notation 7r Ml/ is the space-like traceless part of T ^ v , while the isotropic pressure is the ’’traceful” part. The 
latter is usually divided into the hydrostatic pressure V eq ., that is the pressure that the system would have at global 
equilibrium having £ as the energy density, and the bulk viscous correction II. Therefore 

7r ^ = T <M ) n=-^A ^T^-Peq., (12) 

with the angular bracket we understand the contraction with A(“), defined as 

Kp = 2 ( A « A d + A « A £ - > (13) 

that is, projecting two indices on the spatial, symmetric and traceless part. 

It must be noted that both £ and V eq . are given solely by the local equilibrium background of the distribution, 
while the pressure corrections and II depend only on 5f , 

n = J dP{p ■ A • p)6f , 

therefore, the local equilibrium distribution does not contain any information on the pressure corrections, even if it 
reproduces the full proper energy density. 

More generally, the decomposition (11) is an organization of the ten independent degrees of freedom of the symmetric 
tensor T^", namely the proper energy density £ (and hence the hydrostatic pressure V eq ., depending only on £), three 
independent components of the time-like vector five independent components of it and the bulk pressure II. 
Since this decomposition is completely general, it can be used for the anisotropic background too. 


dPp^p u) Sf, 


(14) 


B. Anisotropic expansion 

The background of the hydrodynamics expansion does not need to be the local equilibrium. In anisotropic hydro¬ 
dynamics one expands the distribution function around an anisotropic background embedding (at least partially) the 
pressure corrections characterizing the full system 

f(x,p) = f a {x,p) + 5f. (15) 

In the previous works [43, 44, 49] it was assumed that the pressure corrections are not fully described by the anisotropic 
background / a , and in general there should be some (expected to be small) contributions from 5f. However, in the 
present work we propose a different approach, extending the Landau matching to the anisotropic case. 

Because of the success of previous applications of anisotropic hydrodynamics, let us consider, for instance, the 
distribution 


fa{x,p) 


k exp 


A(aO 


p^ v {x)p v 


(16) 


which is a generalization of the Romatschke-Strickland form [53]. It must be noted that the form of the anisotropic 
background (16) is not the only one possible. The extended matching procedure we advocate in this Section can be 
applied to different backgrounds. However, in the present work we consider only this background, assuming that the 
generalized Romatschke-Strickland form provides a good enough approximation of the full distribution function and 
deviations from (16) are small perturbations even if there are large gradients in the system. 

The most general decomposition of the symmetric tensor S**", compatible with the Landau prescription of the 
four-velocity, is 


E^ = U^U” + + E^ v) . 


(17) 
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It is convenient now to use the definition of the projector A A “ / in (8) to write 


'ZZf 11 ' _ 


(tu-ta) U»U V + <j> s gr +^ v) . 


(18) 


In order to always have a real valued f a , taking into account that particles are on-shell, both 0y and 4>u — 4> s must 
be positive. Therefore, in a completely general way, it is possible to rescale the quantities and write 


f a {x,p) = A; exp 


( U»U V + 0 gv v + Pv 


= fcexp 


~v ,m 




(19) 


where 


1 = , = 0 S fliV = 

A 2 A 2 0y - 0s ’ 0c/ - 0s' 


( 20 ) 


In Refs [47, 49] it was preferred to have the 0 scalar multiplying the projector This can be easily done here 

by rescaling the quantities in (16) around (j>u instead of 0c/ — 0s- In this paper the prescription (19) is adopted for 
mathematical convenience, the reason will become apparent later. We note that we have now a very similar situation 
to the previous case. Indeed, we recover the usual local equilibrium distribution (5) in the limit of vanishing anisotropy 
tensor E l ,lv . vanishing 0 and, therefore, A becomes the effective temperature T. 

Having fixed the form of the background distribution, we should clearly specify the independent degrees of freedom. 
For the local equilibrium distribution (5), it was necessary to define four parameters, the effective temperature T 
and three independent components of the four velocity field t/ M . Now we aim to generalize the procedure to the 
anisotropic background (19). The four-velocity can be defined adopting the Landau prescription (6), exactly as in the 
previous case. Therefore Sf must fulfill a condition which is equivalent to the one seen for the viscous expansion (7). 
Unfortunately there is no direct way to fix each parameter separately. However, looking at the Landau matching (9) 
for the definition of the temperature, the most straightforward extension is to fix the remaining parameters in such a 
way that the leading order integrals match with the exact values of the whole stress-energy tensor. 

The Landau prescription takes care of the three degrees of freedom for the four-velocity, hence, there are seven 
remaining independent degrees of freedom, namely the momentum scale A, the bulk parameter 0 3 and five independent 
components of the anisotropy tensor . This is the same number of effective degrees of freedom as in the general 
decomposition (11), where we have: the four velocity U 1 0 the proper energy density E , the isotropic pressure V eq . +n 
(and hence n, since V eq . is already defined by E if there are no additional conserved charges), and n^ v . Therefore, as 
a generalization of the Landau matching to the anisotropic case, we are assuming that all of these quantities will be 
reproduced by the anisotropic background alone 


T^ = J dP V > V f = jdP j/Y fa 


( 21 ) 


or, equivalently 


JdPr P v 6f = 0 , ( 22 ) 

which is a generalization of Eqs. (7) and (10) to the anisotropic background. The reason of the choice of the 
parametrization (19) is clear now. In the massless case the system is exactly conformal and the energy-momentum 
tensor is traceless, = 0. Thus, there is one less degree of freedom in . On the other hand, in the anisotropic 
distribution (19), the bulk degree of freedom 0 vanishes too because of the m 2 factor multiplying it. Hence, we still 
have the same number of degrees of freedom in the anisotropic background as in the exact energy-momentum tensor. 


III. DYNAMICAL EQUATIONS 

In this Section we will introduce the equations of motions for the anisotropic expansion. We follow mainly the 
approach of Refs. [25, 26], generalizing the arguments to the anisotropic expansion. The main idea is to get the 


3 It will be shown later that the bulk viscosity IT is proportional to 0 at first order in small deviations from local equilibrium. 
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equations directly from the exact kinematic evolution, then making an approximation on Sf (or its moments) in order 
to have a finite number of equations. 

In order to have energy and momentum conservation, it is necessary to take the first moment of the Boltzmann 
equation (3). Making use of (11), and taking the projection parallel and orthogonal to the four-velocity, one can write 

DE = - ( £ + V eq . + n) 0 + (23) 

(s + v eq . + n) Du a = -v a ( v eq . + n) - a(24) 

where we used the the notation for the convective derivative D = and for the spatial gradient V p = A^dv The 

scalar expansion 0 = dpU p and the shear stress a = A 7 d a U 13 are the ones appearing in the familiar decomposition 
of the expansion tensor 0 pl ' 


0^ = V P U V = cr^ + 0A^ 


(25) 


with = V^t/" — \7 V U P being the vorticity. 

Containing only one convective derivative, Eqs. (23) and (24) are usually considered the time evolution equations for 
the proper energy density £ and the four-velocity U p . It is possible to obtain from kinetic theory a set of equations for 
the remaining degrees of freedom, however, in general it is not guaranteed that the convective derivative will appear 
only acting on n fl1 ' and II. Using the orthogonality properties of the projectors A^, A p and the Landau prescription 
for the four-velocity (6), the exact convective derivatives of the pressure corrections read 

= J dPp V> Df, DM = -1 J dP ( p ■ A • p) Df - DV eq .. (26) 

It is interesting to note at this point that the four-momentum conservation equations and the last equations put 
together are just another way to rewrite the exact evolution equation of all the components of the energy-momentum 
tensor 




dPpPp v Df. 


(27) 


Indeed, making use of the general decomposition (11) and the local four-momentum conservation constraint of the 
collisional kernel J dPp p C[f) = 0, it is straightforward to prove that the UpU v contraction correspond to Eq. (23), 
the A“[/„ one corresponds to Eq. (24), and the remaining independent ones, namely the spatial traceless part and 
the spatial trace of Eq. (27), correspond respectively to Eqs. (26) 4 . 

It is convenient now to write the evolution of pressure corrections in terms of moments of the distribution function 
instead of moments of the derivatives of the distribution function. From the Boltzmann equation (1) it is possible to 
obtain the exact convective derivative of the distribution function 


Df = 


1 

Wu) 


(p-U)U^f = p-df-p-Vf 


C[f \ V ■ V/ 

(p -U) ( P-U )' 


(28) 


Making use of the notation 


= J dP{p-U) r p^...p^C\f}, 


the dynamical equations (26) read 


+ ci7 } = -A£V a J dP pP ^ p ^ f - ( + \e A pff ) J dP 


DU — — A (UI/C7 — DP eq + -A p „Vp 


dP 


p^pVpP f 

( P-U) 


1 


expo T ^ 0 Ap a 


dP 


p^p v )pPp a f 

(p ■ U) 2 ’ 

(p ■ A ■ p)p p p a f 

(p • uy 


(29) 

(30) 

(31) 


Note that no expansion of the distribution function has been used so far, the arguments are completely general. 


4 
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It is important to note at this point that Eqs. (30) and (31) do not depend only on the background f a . In fact, even 
if the second moment of the full distribution / does not depend on the correction Sf because of (22), the same cannot 
be said of the other moments of the distribution (including possible powers of (p ■ U)). Indeed, both the integrals of 
the collisional term on the left-hand side of (30) and (31), and the integrals on the right-hand side cannot be reduced 
to a function of the elements of T in the most general case. Equations (23), (24), (30) and (31) can form a close 
system of equations only at the leading order of the anisotropic expansion, that is, neglecting the contribution of Sf. 

In order to have a more accurate description, it is possible to use the same procedure we used for obtaining Eq. (27) 
(namely, four-momentum conservation plus Eqs. (30) and (31)) to all of the other moments of the distribution function 


D j dP(p-U) r p^ ...p"" f = r DUv J dPip-U) 1 - 1 ^ 1 ...p^p v f - V„ j dP (p ■ U) r ~ x p^ ... p^ n f 

+ (r -1) ( a pa + l -e A pa ^j jdP(p • uy~ 2 yw../y//. ( 32 ) 

This unfortunately provides an infinite set of coupled equations, since the equation for any moment contains higher 
moments of the distribution function. This is not unexpected, indeed, there is exactly the same problem in the 
ordinary viscous expansion [25] around the local equilibrium, once we take into account the different notation. The 
usual procedure is to express the deviation from the background Sf explicitly as series of irreducible moments of Sf 
itself, truncating the series in order to have a finite number of degrees of freedom. See for instance Ref. [50], for 
an extended discussion of different truncation schemes, including different choices of dynamical equations instead of 
the exact moments (32) for closing the set of equations. The expansion in Ref. [50] is around local equilibrium, but 
the arguments can be straightforwardly extended to the anisotropic background. In the specific case of anisotropic 
hydrodynamics, Bazow at al. [43, 44] have been using an “Israel-Stewart like” expansion for Sf. For the reminder of 
this work we will not address the most convenient truncation for the deviation from the background, leaving it for 
further research. Starting from the next Section we will consider only the leading order of the anisotropic expansion, 
namely /„. 

As a final note, as long as we consider the leading order of anisotropic hydrodynamics, the condition (22) is not 
essential, since Sf is systematically neglected. We could still use Eqs. (23), (24), (30) and (31) for closing the system 
of equations. However, in order for the next to leading order to stem from the exact evolution of the moments 
of the full distribution (32) like the leading order, it is mandatory. In fact in Refs. [43, 44, 49], only the Landau 
prescription (7) and the Landau matching (10) were considered, instead of the full generalized matching (22). The 
unfortunate consequence is that Eqs. (30) and (31) can not be used to close the system of equations both for the 
degrees of freedom of the anisotropic distribution function and the extra corrections from Sf . Lacking any clear 
prescription to find the extra needed equations, it comes back the problem of ambiguity in the derivation of fluid 
dynamics already seen in viscous hydrodynamics (see the discussion on Israel-Stewart approach in Ref. [50]). For 
instance, in Ref. [43, 44] the needed extra equations for closing viscous anisotropic hydrodynamics have been selected 
from the zeroth and second moment of the Boltzmann equation. 


IV. SPECIFIC EXAMPLE: LEADING ORDER FOR THE RELAXATION-TIME APPROXIMATION 

In this Section we discuss the specific case of the collisional kernel C [f] treated in the relaxation-time approximation 
(RTA) 


C[f] = -( p-U ) f(x,p) - f oq .(x,p) 


(33) 


with r eq being the relaxation time and / eq . the local equilibrium distribution (5) 5 . 


5 The local equilibrium distribution is well defined in any situation through the four velocity prescription (6) and the Landau matching (9), 
even if it is not used as the background in the hydrodynamics expansion. 





Thanks to Eq. (33) the collisional kernel in Eqs. (30) and (31) can be easily calculated, and we can write 


+ = “ A - v " / ipt ww - (■ 


r pa + 3 S A p<j 


dP 


(n + p eq .) 


i 


l 


D ( n + Vea. ) h -n — + — A uv v 


1 eq. 


, i , 


per p;- ,, ^ A po . 


dP 


p(»p-')pPp° f 

(p • CO 2 ’ 

(p ■ A ■ p)p p p CT / 

(P • c/) 2 


(34) 

(35) 


We are interested in the leading order of the expansion, therefore we will consider the approximation / ~ f a with the 
anisotropic distribution function defined in (19). Because of the special form of the collisional kernel in (33) and the 
extended matching (22), the left-hand sides of both Eq. (34) and Eq. (35) do not depend on Sf. On the other hand, 
the first term on the right-hand sides of both equations can be written directly in terms of the components of T pL/ . 
Indeed, because of the symmetry of f a 6 


dP pPp P r J a = jj p T aa + U a T ap + U a T pa -2 £ U p U a U a 


(p-U) 


therefore, Eqs. (34) and (35) read 


Dtt^ + —7 = - 
Teq. 


£>( n + p eq .) 


Tp* + \o&p^j J 
+2 ( P eq . + n) - ^9 + 2 7 T<»u v>a , 

(P ' A • p)p P p a f a 


+ T ca n ~ 3 ( apa + 3 dApa 


dP 


C p-u) 2 


(36) 


(37) 


(38) 


+ 2 — g (’Peq. + II) 9 . 


It is convenient at this point to introduce the notation for the isotropic pressure 

1 


v = v cq . +n = - 3 


J dP (p-A-p)f a , 


(39) 


since, differently from the hydrostatic pressure P e q. itself, it does not depend on the equilibrium distribution function 
/ eq .. Equations (37) and (38) then read 


Dn^ + —-— 7r pL/ = - 


j + - 9 A ) [ dP ^^ fa - 2 ^a v > a 

pa + 3 ° pa J J a (p-U )2 Z7r a a 


+2 V a pv - -6 + 2 7 T< p u v>a , 

o 


DV 


+ Teq ( v Peq.) - 3 ( + 3 6 A P-) / 


dP 


(P ' A • p)p P p a fg 

(p • uy 


(40) 


(41) 


+\w ,tv -\vo- 


A. Close-to-equilibrium limit 


It is well known that second order viscous hydrodynamics is a very good approximation if the system is close to 
local equilibrium. Therefore, if the prescription we present for anisotropic hydrodynamics is a good description of 
the system, we expect that our equations should match with the second order viscous hydrodynamics in the close-to- 


6 In fact if Sf preserves the parity invariance (seen from the local rest frame) of f a , this is still valid. However in the most general case 
there are some other non vanishing terms. 
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equilibrium limit. Namely, Eqs. (37) and (38) should reduce to 7 

tvDtt^ + tt^ = 2r] a 111 ' + 2 r n 7t^ m a v ^ x — 6^ + A^n n a y,v , 

m MI + II = — (d — ^nn 110 + An^ 


(42) 

(43) 


with certain values of the transport coefficients. 

In the previous works on anisotropic hydrodynamics it has been proven that the correct close-to-equilibrium be¬ 
haviour is reproduced, see for instance Refs. [31, 45, 49]. To demonstrate this property, it was necessary to make 
an expansion of the anisotropic background around the local equilibrium, keeping only the lowest order terms, and 
checking that the equations obtained in this approximation reproduce the second order viscous hydrodynamics. In 
the present work such checking of the close-to-equilibrium limit is easier, since the equation of motion for the pressure 
corrections (34) and (35) are exactly the same as the two prescriptions for second order viscous hydrodynamics tested 
in Ref. [52]. As we show below, the linear expansion of f a around equilibrium yields exactly the distribution function 
known from the gradient expansion. Therefore, the transport coefficients in (42) and (43) are exactly the same as 
those obtained in the Chapman-Enskog-like expansion presented in Ref. [52], which is in our opinion the best second 
order viscous hydrodynamics prescription for dealing with the RTA kinetic theory. 

Indeed, if we take the expansion of (19) up to the leading order in the deviation from local equilibrium we have 


fa 


/eq. 


1 + ±-^(p.U) 


m 2 <j) 

2T(p ■ U ) 


1 

2 T(;p • U) 


P-£-P 


(44) 


Introducing the notation 


In ’ q ~ {2q + l) 


jj/ dP (p • U) n ~ 2q (—p • A • p) q f oq ., 


(45) 


the energy-momentum tensor then reads 




A-T T m \ T 

rp 2 1 3,0 - ^ V 1 1,0 


U^U V - 


A - T T m 2 

p 2 7 3,i ~ 2 T 




(46) 


where X , 'i" stands for 


T£ = J dP P y'f eq . = £U»U V - V e q.A^. (47) 

Since the energy density cannot have any corrections, because of the definitions of both / eq . and / a , we can remove 
A from Eq. (46). Indeed, since 


A-T T m \ T 

1 3,o — Tpyr <p -*i,o 


2^2 


2 T 


A -T _ mV Jqo 
T 2 ~ 2T 0 J 3j0 ’ 


the energy-momentum tensor reads 


rrifll/ ^ rrifliy _ 

- eq. ^ 


-*1,0 

L-^3,0 


I 3.1 ~ h. 




We can recognize then 


= -7 F h,2e v = -P«? iv , 


T 
2 

„ m 
= 2T ' 


1 ’° T T 

- - ^3,1 — -*1,1 


L^3,0 


(48) 


(49) 


(50) 

(51) 


7 Note that in general there can be more second order terms, we are only showing the non vanishing ones in the close-to-equilibrium limit 
of anisotropic hydrodynamics. 
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where is the same as the coefficient in Ref. [52]. Making use of the last equations and of Eq. (48), expansion (44) 
reads 


fa - /eq. < 1 - 


n 


{p-U) 


h,0 T T 


1 -^(P-U) 2 

I 3,0 


2 Th 


P-TT-P 


Using the exact relations between the I n ^ q thermodynamical integrals 

1 


In. a — 


2q + ll 


In,q— 1 VII In—2,q—l 


T =T 

*n,q 1 


n—l,q— 1 


( n - 2q^j 


n—1 ,q 


the expansion of the anisotropic bakground (52) reads 

n 


f a ~ /eq. 1 - 


3 T(p-U)k 


h,i 


m z - 1-3-^ (p-U) 


*3,0 


+ 


2 T/3, 


- p ■ n ■ p > , 


where 


Pn = ~ ( ^2,0 + ^2,1) ■ 

•j *3,0 v / 


(52) 

(53) 

(54) 

(55) 

(56) 


The right-hand side of (55) matches exactly the form of the distribution function 
and (35) for obtaining the second order viscous hydrodynamics in Ref. [52], q.e.d. 
possible to plug directly the expansion (44) into the right-hand side 8 of Eqs. (40) 
results. 


that was plugged in Eqs. (34) 
Naturally, it would have been 
and (41) and obtain the same 


B. Numerical results 

In this Section we present direct comparisons of the new prescription for the anisotropic hydrodynamics expansion 
with the exact solutions of the Boltzmann equation. We consider the exact solution found in Ref. [42], namely the 
0 + 1-dinrensional solution for a massive gas with the collisional kernel treated in the relaxation-time approximation. 
We consider a fixed relaxation time r eq = 0.5 fm/c, a starting temperature T 0 = 600 MeV, and a starting time r 0 = 0.5 
fm/c. For the mass we consider two options, namely, m = 0.3 GeV and m = 1 GeV. For the initial anisotropy, we take 
into account two extreme situations, the system that is initially in local equilibrium, £o = 0, and the system with a 
very large starting anisotropy, = 100. In Appendix VI we present the 0 + 1-dinrensional expressions of the equations 
presented in this paper. We plot the evolution of the pressure anisotropy, that is the ratio of the longitudinal pressure 
Vl and the transverse pressure 'Pt 9 , and the bulk pressure correction times the proper time rll. 

We show the comparison with the exact solutions of three different approximations. By anisotropic matching we 
understand the leading order of the prescription presented in this paper, by viscous hydrodynamics we denote the 
modified gradient expansion by Jaiswal et al. [26] and, finally, anisotropic hydrodynamics stands for the prescription 
introduced by Nopoush et al. in Ref. [47]. The latter assumes the longitudinal boost invariance and cylindrically 
symmetric radial expansion. There is a recent, more general, formulation of anisotropic hydrodynamics [49] which 
does not require such symmetries in the evolution of the system. However, as it is shown in Ref. [65], the less general 
prescription has a similar or even better agreement than the more general one in the (0 + l)-dimensional framework. 
On the other hand, it is useful to remind that the viscous hydrodynamics prescription we choose, is the close-to- 
equilibrium limit of the equations seen in the last Section, therefore, we should expect a better agreement between the 
full equations and the exact solutions of the Boltzmann equation, than between the linearized version and the exact 
results. Confronting the last results found in the literature (see, in particular, Refs [51, 52]), it is interesting to note 
that the prescription by Nopuosh et al. for anisotropic hydrodynamics and Jaiswal’s prescription for the second order 


8 The left hand side is already exact at any order in the expansion of the distribution function. 

9 See the Appendix for more details. 
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fo=0, m=0.3 GeV 



fo=100, m=0.3 GeV 



r [fm/c] 


fo=0, m=0.3 GeV 



fo=100, m=0.3 GeV 



r [fm/c] 


FIG. 1: (Color online) Time evolution of the pressure anisotropy Vl/Vt (left) and the bulk corrections time the proper time 
rll (right). The bleck solid line correspond to the exact solution of the Boltzmann equation, the red dotted one to the set of 
equation presented in this paper, the purple dashed line is the solution of second order viscous hydrodynamics, while the blue 
dashed-dotted is anisotropic hydrodynamics. The initial anisotropy is vanishing (top) or £o = 100 (bottom) in the two rows. 


viscous hydrodynamics provide the best agreement with the exact solutions, excluding higher order treatment 10 . In 
particular the anisotropic hydrodynamics prescription systematically reproduces the pressure anisotropy better, while 
the viscous hydrodynamics one is more accurate for the isotropic pressure corrections. 

Figures 1 and 2 clearly demonstrate that the approach providing the best approximation to the exact solution is 
the leading order of the anisotropic matching procedure, i.e., Eqs. (23), (24), (40) and (41) applied to the (0 + 1)- 
dimensional expansion. The anisotropic matching provides a similar or better agreement for Vl/Vt than anisotropic 
hydrodynamics, and at the same time a similar or better agreement for rll than viscous hydrodynamics. In particular, 
for large starting anisotropy, it removes the qualitatively wrong behavior of the bulk evolution at initial stages, which 
has been up to now a common feature of leading order anisotropic hydrodynamics, see for instance Ref. [65]. At last, 
as it is expected, the equation we propose reproduce the exact results better the viscous hydrodynamics prescription 
(which we remind correspond to the close to equilibrium limit of the anisotropic matching), especially for the pressure 
anisotropy. This suggests that the higher order corrections are properly taken into account using the new procedure 
already at the leading order. 


10 For the next to leading order treatment of anisotropic hydrodynamics see Refs. [43, 44], and for viscous hydrodynamics see Ref. [43, 66]. 
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£o=0 , m=1 GeV 



£o=100 , m=1 GeV 



£o=0 , m=1 GeV 



fo=100 , m=1 GeV 



FIG. 2: (Color online) Same as Fig. 1, except here we take m = 1 GeV. 


A very interesting comparison can be made between Fig. 2 of the present paper and the plots having the same 
parameters in Ref. [44] 11 . The agreement with the exact solutions for the anisotropic expansion (still at leading order 
in the present paper) is already very close to that found for the viscous anisotropic hydrodynamics , that is, next 
to leading order anisotropic hydrodynamics, but without making use of the generalized Landau matching (22). It 
would be interesting to check in the future if this striking agreement with the exact solution is preserved in other 
dynamical situations, for instance, for different masses or different flow configurations ( e.g. anisotropic hydrodynamics 
has already been tested for the Gubser flow [48]). 


V. SUMMARY AND CONCLUSIONS 

In this paper a novel method to treat the hydrodynamics expansion around an anisotropic background is intro¬ 
duced. Differently from previous works on this topic, we propose an extended Landau matching for the anisotropic 
distribution. In other words, we assume that the anisotropic background is capable to reproduce the stress energy- 
momentum tensor, without the contribution of the higher order corrections (21). In the present work only the 
generalized Romatsclrke-Strickland form (19) is considered as the leading order distribution, but our arguments can 


11 That is Fig. 1. Note, however, the log-log scale in 'Pl/'Pti ami the different prescription for viscous hydrodynamics. 
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be easily adjusted to a different anisotropic distribution, if needed. 

The use of the generalized matching allows us to take the dynamical equations of hydrodynamics directly from the 
exact time evolution of the moments of the Boltzmann distribution (not to confuse with the moments of the Boltzmann 
equation). Following this procedure it is necessary to include at least the evolution of the energy momentum tensor 
(second moment of the distribution), since Eq. (27) includes local four-momentum conservation (3). Even if the latter 
equations depend only on the stress-energy tensor itself (and, thus, on the leading order only, because of the generalized 
Landau matching), the other independent equations get some contribution from the full distribution function because 
of higher-orders moments. This is a common feature of the exact evolution of all moments of the distribution function, 
namely, some higher moments of the distribution always appear in their exact evolution equations (32). Therefore, 
the exact dynamics is given by an infinite set of coupled equations for the moments of the Boltzmann distribution. 
The solution of this puzzle is equivalent to the one used in viscous hydrodynamics, that is, one should make an 
approximation for the deviation from the anisotropic background having a finite number of degrees of freedom (which 
are connected directly with the moments of the full distribution itself). In this way one can, for instance, continue 
using the exact evolution of the moments of the distributions to close the system. 

In the present paper we consider for practical applications only the simplest approximation, restricting ourselves 
to the leading order. Namely, the full distribution is assumed to be well reproduced by the anisotropic background 
without the need of higher order corrections. As a consequence, Eq. (27) is enough to close the system of equations. 
We leave the selection of the most appropriate treatment of the next to leading order corrections to be addressed in 
further research. 

We analyze in detail the case of the collisional term treated in relaxation-time approximation and prove that in the 
close-to-equilibrium the equations of motion reduce to the most accurate prescription for the second order viscous 
hydrodynamics. We find a striking agreement between the new approach and the exact solution of the Boltzmann 
equation in the Bjorken flow limit. Both the pressure anisotropy and the bulk evolution are significantly improved 
compared to the most successful approximation; namely Nopoush et al. version of anisotropic hydrodynamics for 
the pressure anisotropy, and Jaiswal modified Chapman-Enskog viscous hydrodynamics for the bulk [51, 52]. These 
preliminary results suggest that it is necessary to take higher orders approximations to find a better agreement with 
the exact solutions. For instance, next to leading order anisotropic hydrodynamics or third order viscous anisotropic 
hydrodynamics, as the ones in Refs. [43, 44]. In fact, comparing the plots having the same parameters in the present 
work and in Ref. [44], the agreement with the exact solutions already seems very close to next to leading order 
anisotropic hydrodynamics. 
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VI. APPENDIX: EXPLICIT FORMULAS 

Here we present the formulas used for the numeric calculations in the last section. Because of longitudinal boost 
invariance and homogeneity on the transverse plane, the four velocity reads for the Bjorken flow 


(57) 

being ?y the space-time rapidity. All Lorentz scalars depend only on the longitudinal proper time 


U M = (cosh 77, 0 , 0 , sinh ?y), 


ry = -ln 


t + z 
t — z 


T = \/£ 2 — Z 2 . 

The scalar expansion and the shear stress therefore read 



2 

3 r 




(58) 


^ ^ ( X»X V + Y^Y V 


(59) 
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where 


Z^ = (sinhjj,0,0, coshr/), (60) 

is the longitudinal direction, and X ,J ’ and Y^ are two arbitrary directions on the transverse plane, for instance the x 
and y direction in the lab frame. The vector fields U 11 , X and Z M provide an orthonormal basis, which in the 
local rest frame reads 


£7" = (1,0,0,0), ^ = (0,1,0,0), Y* 1 = (0,0,1,0), Z" = (0,0,0,1). (61) 

Because of the symmetry, even the distribution function (19) is largely simplified in the 0 + 1-dimensional setup. 
Indeed, seen from the local rest frame (LRF) 


fa = k exp 



(62) 


being p l = (p • Z) the longitudinal momentum, and pt the transverse momentum. It is convenient at this point to 
rewrite the last formula in the following way 


fa = k exp 



(63) 


where 


1 ! + (/> 1 _ 1 ~ K 1 _ l + £ 

A 2 A 2 ’ (1 + </>)’ af (1 + cf>) 


(64) 


Just like for the shear stress er^", the symmetry of the expansion constrains 7to depend on a single scalar function 
7T S 


TT^" = ^7T s (r) ( X»X V + - 7 T a {r)Z»Z v . (65) 

It is convenient to introduce the the longitudinal pressure Vl and the transverse pressure Vt> 

V L = Z-T-Z = V eq .+H-n s , Vt = X-T-X = Y-T-Y = V eq .+n+K s , (66) 

since this results in relatively simple expressions, which do not depend on the effective temperature, namely 


V L = 

k 1 

f W3p PL 

J 

yjm 2 -f p 2 

Vt = 

k 

[j 3 r Pt 

2 j 

1 \/m 2 - l-p 2 


expf-I, L 2 + 4 + 4)=^A 4 , 


A' 

1 


X T 


e x p|- T1 /m^+?| + 5|)=M, T A 4 . 


A 


(67) 

( 68 ) 


Here the integrals are taken in the local rest frame. Having the same anisotropic distribution for the same (subset of) 
equations, the energy and momentum conservation reduce to the single equation originally found in Ref. [47], namely 


DS = — 


£ + V L 


(69) 
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or, writing explicitly the derivatives on the variables A, cut and 


( 4H 3 - f\,.) 


4d~Lz — n m ) IllnA + 2VLt Dinar + —h D In ax = 0. 


1 


(70) 


The explicit expressions of the functions 'Hz, Hzt, Hzl, and f 1l are the ones found in Ref. [47] itself. They 

are all derived from the exact formula for the proper energy density in 0 + 1-dimensions 


£ = A 4r Hz = {2n k ) a? A 4 / dr r 3 exp 

Jo 


-St 


H 2 (y,z), 


being 


m cxl 

m = T , y= —, 

A OtT 


z = 


m 

arr ’ 


and 


H 2 (y,z) = y J d{— cos 6)^'z 2 + y 2 cos 2 9 + sin 2 6 = 


Vv 2 ~ 1 


(1 + z 2 ) tanh 1 \{ y . ? , \ +\/ ( V 2 + z 2 )(y 2 - 1) 


y z + z 


The other other functions are 


Hz l = (27r k)a^ / dr r 3 exp -\Jvn? + r 2 H 2 L (y,z), 

Jo 1 J 

Y /-OO _. 

Hzt = -(27t fc)ttT / dr r 3 exp — \/di 2 + r 2 
2 Jo 

Hzm 


H 2t {jj, z') , 


(2n k)a^m 2 / dr 

Jo 


y/m 2 + r 2 


exp 


— \J to 2 + r 2 


H 2 {y,z), 


while 


H 2L {y,z) = yd y U 2 {y,z) - H 2 {y,z), 

H 2T (y,z) = H 2 (y, z) — 'H 2L {y, z) — z8 z H 2 (y, z), 


and 


(71) 


(72) 


(73) 


(74) 

(75) 

(76) 


(77) 

(78) 


BH 3 

Bar 

dUz 

Bolt 

dHz 

Bolt 


2 

CUT L 
1 

a L l 
1 ' 
m . 


= —^r, 

OLT 


'Hz + 'Hzt 
'Hz + Hz l 
H z — Hzl — 2.H.ZT ~ H 3 n 


= - &L, 

aL 


m 


(79) 

(80) 
(81) 


Taking into account the symmetry of (63) and making use of (59), (65) and (66), we can rewrite equations (40) 
and (41) in the following way 
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D\ 


1 (Vt ~ Vl) = -A 4 ^ H.3TL — P^Ll) - 
( 2 Vt — Vl) = -A 4 ^ 2ii 3 TL + P-3 Ll) 
The functions H 3 tl and H 3 ll are defined by 


D 


- f- + —) ( 2 V T + Vl) - ~V L + —V, 

\ T ^"eq. / ' / T 7~eq. 


(82) 

(83) 


P-3LL 

P-3TL 


1 

A 4 

1 


W dP 


dP Jp^jf fa = ^ 27Tk)a ^ 

(l vX)\p-Z)\ _ 1 


dr r 3 exp 


— \/ m 2 + r 2 


P 


2LL 5 


(p-t /) 2 


f a = ~(2irk)a T 


dr r 3 exp 


-v^ 


to - 


P 2 TL, 


being 


(84) 

(85) 


'H2TL = %2T — ydy'H. 2 T, P 2 LL = 31-L2L — ydy'H. 2 L- (86) 

In order to solve Eqs. (82) and (83) it is necessary to have an explicit formula for the derivatives of the longitudinal 
and transverse pressure 


being 


DV l = A 4 
DV t = A 4 


( 4P-3L — &Lm) D In A + VLltD In olt + &llD In , 
^ / PH. 3 T — &Tmj /I hi A — &ttD In qj- + CItlD In , 


^2 Lm — 7 / 7 . dm P 3 l ) 

^Tm — 777 dm P 3r , 


ilLT = OtT d aT P 3l , 
&TT = OtT d aT P 3t , 


Introducing the quantities 


Q.LL = OIL da L 'H 3l , 
VLtL = OiL da L H 3t . 


P-3Lm = (27T k) 4 / dP 


y/rh 2 +: 


: exp 


— y/ rh 2 + r 2 


U 3T m = -(2tt k)a% / dP 


m 


y/rh 2 + \ 


: exp 


— y/ rh 2 + r 2 


P-2L, 

P 2 T, 


Pstt = -(27rlc)4 f dP exp 


— yj rh 2 + r 2 


P-2 TT, 


with 


(87) 

( 88 ) 


(89) 

(90) 


(91) 

(92) 

(93) 


W-2TT = ydy P 2 T + z d z V.2T, 


(94) 


it is possible to write all of the fl’s functions as integrals of the derivatives of (73). 
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&Lm = 'H 3 LL + ‘2'HzTL ~ L ~ (95) 

VLll = 3 H 3 l — T-L^ll, (96) 

VLlt = ZHsl — 2 H 3 tl, (97) 

9rm = 27^3 TT + 773TL — 'H 3 T — 7^3 Tm, (98) 

(7tL = 773T — i~i-3TL, (99) 

( 2 tt = 4773T — 2H3TT- (160) 
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